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Efimov effect from functional renormalization 
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Institut fiir Theoretische Physik, Philosophenweg 16, D-69120 Heidelberg, Germany 

We apply a field-theoretic functional renormalization group technique to the few-body (vacuum) 
physics of non-relativistic atoms near a Feshbach resonance. Three systems are considered: one- 
component bosons with a U{1) symmetry, two-component fermions with a U(l) x SU{2) symmetry 
and three-component fermions with a ?7(1) x SU{?>) symmetry. We focus on the scale invariant 
unitarity limit of an infinite scattering length. The exact solution for the two-body sector is consis- 
tent with the unitary fixed point behavior of the considered systems. Nevertheless, the numerical 
three-body solution in the s-wave sector develops a limit cycle scaling in case of U{1) bosons and 
SU{2,) fermions. The Efimov parameter for the one-component bosons and the three-component 
fermions is found to be s « 1.006, consistent with the result of Efimov. 

PACS numbers: 21.45.-v; 34.50.-s 



I. INTRODUCTION 



The physics of ultracold atoms is a broad area of research which develops rapidly both experimentally and theoreti- 
cally (for reviews see [ll,[2|)- To a large extend this is due to the excellent tunability and control of the studied systems, 
"j^ \ In particular the interaction strength of atoms near a Feshbach resonance can be changed in broad ranges by tuning 
the magnetic field, which makes these systems an ideal playground for testing the predictions of theoretical models at 
strong coupling. Both, few-body and many-body quantum, and thermodynamic effects have been extensively studied 
"j^ , with ultracold gases. 

' Near a broad Feshbach resonance the atom-atom s-wave scattering length a is large in comparison to the effective 
I , range reS of the microscopic interaction. The low-energy vacuum physics (for vanishing temperature and density) 

■ becomes universal: some physical observables become insensitive to the detailed form of the microscopic interaction 
Ch ' and depend only on the scattering length a For example, for a > the theory admits a stable shallow diatom. 

For this atom-atom bound state the universal binding energy is determined simply by dimensional analysis. In the 

^ unitarity limit all energy scales drop out of the problem and the theory is scale invariant in the two-body sector. It is a 

well-established result, derived first by Efimov [4l|, that in the three-body sector of the resonantly interacting particles 
' a spectrum of shallow three-body bound states develops. At unitarity, the spectrum is geometric which is a signature 
, of the limit cycle behavior of the renormalization group flow. Even in the case of a scale symmetry in the two-body 

■ sector, the running of the renormalized three-body couplings indicates a violation of the dilatation symmetry and 
(y-j \ nray be associated with a quantum anomaly [H] ■ 

■ The low-energy few-body scattering of atoms has been investigated using various computational non-perturbative 
techniques ranging from effective field theory 0, d, Q to quantum mechanics [1, |T0|] . The perturbative e expansion 
around critical d — A and d = 2 dimensions has also been applied to this problem [ll|, Tl^ . A field theoretical 

' functional renormalization group approach has been used to investigate the two-body and three-body sectors of 
two-component fermions recently |13i , ,14] . As a convenient truncation in vacuum, the authors use a vertex expansion 
^ and reproduce the Skorniakov and Ter-Martirosian integral equation 15]. In this way, the universal ratio of the 
. 1-H atom-diatom to the atom-atom scattering length is computed. 

5h In this work we follow pj| and consider the few-body physics of non-relativistic atoms near a Feshbach resonance 

. . . 1 which may be described by a simple two-channel model of particles with short-range interactions. We study three 
different systems: bosons with a U{1) symmetry (System I), fermions with a J7(l) x SU{2) symmetry (System II) 
and fermions with a U{1) x SU{3) symmetry (System III). Both, Systems I and II have been well-studied during the 
last decade. The model of SU (3) fermions might be of relevance for three-component mixtures of ^Li atoms near 
the broad Feshbach resonances. The many-body properties of this model have been studied in [l^, [13, 
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Recently, the three-component fermion system has been studied with functional renormalization group methods using 
an approximation including a trion field (21i] . The present work, which is based on a vertex expansion, complements 
and extends the results of [2l| . It underlines the basic finding of the presence of Efimov states for SU (3) fermions 
and estimates the universal Efimov parameter s with a higher precision. 

The structure of the paper is as follows: In Section [TT] we present a field-theoretic renormalization group (RG) 
method, which we use to solve the few-body problem, and introduce the three models we are going to investigate in 
this work. In Section [IIII we investigate the effective action in the vacuum state, i.e. for vanishing temperature T = 
and density n = 0. The following Section IIVI is devoted to the exact solution of the two-body sector for positive 
scattering lengths a > (diatom phase). In Section |V] we turn to the analysis of the three-body sector and derive the 
RG flow equation for the atom-diatom vertex at unitarity. This RG equation is solved analytically employing a simple 
pointlike approximation in Section IVII In Sections IVIII and I Villi we reproduce the Skorniakov and Ter-Martirosian 
integral equation and present a numerical solution of the three-body RG flow equation. We draw our conclusions in 
Section HXl 



II. METHOD AND DEFINITION OF MODELS 

In this work we calculate a scale-dependent effective action functional F^. (for reviews see @, [13]), often called 
average action, flowing action, or running action. This renormalization group method is formulated in Euclidean 
spacetime using the Matsubara formalism. The flowing action Tk includes all fluctuations with momenta q ^ k. In 
the infrared limit fc — > the full effective action F — Ffc^o is obtained. This dependence on the scale k is introduced 

(2) 

by addinga regulator Rk to the inverse propagator Fj, and the flowing action F^ obeys the exact functional flow 
equation [22l |: 

dkVk = isTrafci?fe(Ff +i?fc)-i = isTr9fc ln(Ff +i?fc). (1) 

This functional differential equation for F^ must be supplemented with the initial condition Tk^\ = S, where the 
"classical action" S describes the physics at the microscopic UV scale, fc = A. In Eq. ([!]) STr denotes a supertrace 
which sums over momenta, Matsubara frequencies, internal indices, and fields (taking fermions with a minus sign). 

(2) 

The second functional derivative F^ denotes the full inverse field propagator, which is modified by the presence of 
the IR regulator Rk- As a consequence, the fiuctuations with < fc^ are suppressed and the effective action depends 
on the scale k. The choice of the momentum dependent regulator function Rkiq) introduces a scheme dependence 
which has to disappear for the exact solution for fc — > 0. In the second form of the flow equation ([T|) dk denotes 
a scale derivative, which acts only on the IR regulator Rk- This form is very useful because it can be formulated 
in terms of one-loop Feynman diagrams. The effective action Ffe=o is the generating functional of the IPI vertices, 
which can be easily connected to the different scattering amplitudes in the case of vanishing density (n = 0) and 
vanishing temperature (T = 0). It is also convenient to introduce the RG "time" t = ln(fc/A), which flows in the 
interval t G (— oo, 0). In the following we will use both t and fc. 

In most cases of interest the functional differential equation ([1]) can be solved only approximately. Usually some 
type of expansion of Tk is performed, which is then truncated at finite order leading to a finite system of ordinary 
differential equations. The expansions do not necessarily involve a small parameter (like an interaction coupling 
constant) and they are, in general, of non-perturbative nature. As has already been advocated in Sect. HI we perform 
a systematic vertex expansion of F^: taking the full momentum dependence of the relevant vertex in the three-body 
sector into account. The vertex expansion is an expansion in powers of fields; hence generally: 

OG 

^k = J2 W = rfe(2) + Ffe(3) + Ffe(4) + (2) 

where the index in brackets denotes the number of fields n in the monomial term Ffe(rt). In the second equation 
Ffe(O) and Ffc(l) are missing because we are not interested in the free energy of the vacuum and the term linear in 
the fields is absent by construction. 

In this paper we are interested in the non-relativistic physics of atoms interacting via a Feshbach resonance, which 
can be described by a simple two-channel model. In particular, we consider and compare three different systems: 

• System I: Single bosonic field near a Feshbach resonance 
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Our truncation of the scale-dependent flowing action, written in the Fourier space, is: 
rfc(2) = / i^*(g)(zc^q + q2-^^)V^(g)+ / ip*{Q)P^{Q)^iQ) 



rfe(3) = J / f*{QimQ2MQ3) + 'p{Qi)riQ2)r{Q3) S{Qi^Q2~Q3) 



Q 

h 
2 

rfe(4) = - J X3iQi,Q2,Q3MQi)^{Q2)v*iQ3)riQ4)SiQi + Q2 - Qs - Qi), (3) 

where Q = (ti^,q) and J = ^ (2tt)-^ • field "0 represents an elementary complex bosonic atom, 

Q °° °° 

while V'(Q) is a complex bosonic composite diatom which mediates the Feshbach interaction. At the initial UV 
scale we take A3 = 0. The action for Lp becomes Gaussian, and one may integrate out Lp using its field equation 
ip ^ ipip. As will be demonstrated in Sect. IIVI the Yukawa coupling h is simply related to the width of the 
Feshbach resonance. For k ^ the coupling X3{Qi, Q2, Q3) becomes the IPI vertex which can be connected to 
the atom-diatom scattering amplitude. The system has an obvious U{1) symmetry which reflects the conserved 
number of atoms^. 

• System II: Fermionic doublet near a Feshbach resonance 



2 

Tk{2) = E/ V'*(Q)(*c^q + q'-Mv)V'.(Q) + J v*{Q)P^{QMQ) 



rfc(3) 

r;o(4) 



Q 

v*{Qi)MQ2)MQ3) - v{Qi)ri[Q2W2{Qz)\5{Qi -Q2~ Q3) 



Ql:Q2:Q3 



A3(Qi, Q2, Q3) ^iQi)MQ2)^*{Q3miQi)SiQi + Q2-Q3-Q4 



(4) 



0i,...Q4 



Here, the two species of elementary fermionic atoms '4'i,'4'2 are described by Grassmann- valued flelds, and (p is 
a composite bosonic diatom. At the UV scale one has (p ~ ipi'4'2- This fermionic system has an SU{2) x U{1) 
internal symmetry with ("01, "02) transforming as a doublet and (p as a singlet of the SU{2) flavor subgroup. 
Two-species fermion systems near Feshbach resonances were realized experimentally with ^Li and *'^K atoms 

• System III: Fermionic triplet near a Feshbach resonance 

3 » 3 



rfe(3) 
rfc(4) 



Q 
h 
2 



3 

^ ey-fek(gi)V',(Q2)V'fe(Q3) - v^{Ql)r,{Q2)rk{Q3)\5{Ql - Q2 - Q3) 



3 3 

A3a(Qi, Q2, Q3) ^^iQi)MQ2) Y ^*jiQ3)^jiQd 



3 3 
A36(Ql, Q2, Q3) Y MQ2miQ4) Y V>jiQl)^*jiQ3)] '^(Ql +Q2~Q3-Q4 

t=l j=l 



(5) 



^ In general, the U{1) symmetry can be spontaneously broken due to many-body effects and our truncation Jsjl would be insufficient. In 
this work, however, we are interested only in the few-body physics (for more details see Sect. IIIII l. 



4 



The three species of the elementary Grassmann-valued fermion field can be assembled into a vector ^|J = 
(V'l) V'2j V's)- Similarly the three composite Feshbach bosonic diatoms form the vector tp — {fi,ip2iV3) ~ 
V'aV'ii ■'/'i'02)- The action has an SU{3) x U{1) symmetry with ip transforming as 3, and as 3 for the 
SU{3) flavor subgroup. Two different couplings Asa and Ash are allowed by the SU{3) symmetry. This model 
might be of relevance for three-component mixtures of ^Li atoms. There are three distinct broad Feshbach 
resonances for three scattering channels near B « 800 G for ^Li atoms. As a first approximation wc assume that 
the resonances for all channels are degenerate, which leads to the SU{3) fiavor symmetry and to the model (O. 
A stable three-component mixture of ^Li atoms has been recently created ^26, 27]. The theoretical investigation 
of the 3-body losses in [2^, [23| has been recently published [H, IH, [s^l ■ 

To unify our language for the different models considered in this paper, we refer to the elementary particles tp 
as atoms (and denote corresponding quantities with the subscript ip), while the composite tp is called diatom. All 
considered systems have Galilean spacetime symmetry, which consequences we discuss in Appendix [A] Our units are 
h ~ kB ^ I- Moreover we choose the energy units such that 2M^ = 1, where M^f, is the mass of the atom. 

We should stress that r(2) and r(3) do not have the most general form. The most general form of the ver- 
tex expansion includes an arbitrary inverse atom propagator and a momentum-dependent Yukawa coupling 
(527 Qs)- However, due to special properties of the vacuum state (see Sect. IIII|) . P^{Q) and /i(Qi,Q2jQ3) are 
not renormalized and keep their microscopic values P-ip{Q) = (iti-'q + — Ati/i) and h{Qi, Q2, Q3) — h during the RG 
fiow. At this point it is also important to note that our vertex expansion is complete to the third order in the fields. 
Possible terms with four fields, which are invariant with respect to the symmetries of our models, can be found in 
Appendix |B] In this Appendix we also present arguments, because of which we do not include these terms in our 
truncation. To summarize, the properties of the two and three-body sectors, which are of the main interest in this 
work, can be calculated using the truncations IS]) , ([4]) , ([S]) . 



III. VACUUM LIMIT 



The advantage of the method used in this paper is that it is a field-theoretical setting which permits computations 
for the general case of non-zero temperature (T 7^ 0) and density (n ^ 0). In this work we are interested only in 
the scattering and the bound states of few particles in vacuum. The projection of the effective action rfc=o onto the 
vacuum state must be performed carefully and was developed in l3l|. Here we shortly summarize the procedure: 

The vacuum projection of F^^^o is performed as follows: 



r^ac — lim Tk=o 



T>T,{kF) 



(6) 



where kp = (37r^n)^/"^ is a formal Fermi wave vector (defined for both bosons and fermions) and n is the atom 
density of the system. Thus we start with the effective action at finite density and temperature. The system is then 
made dilute by taking limit kp — > 0. It is crucial, however, to keep the temperature T above its critical value in 
order to avoid many-body effects (e.g. Bose-Einstein condensation). One may perform the vacuum limit for a fixed 
dimensionless ^ such that the temperature goes to zero because Tc scales ~ fc|n. 

Let us now examine the momentum- independent part of the atom inverse propagator P^^k=o{Q = 0) = — as 
well as its diatom counterpart = P^,k=o{Q = 0), in more detail^. For positive values, i. e. fj,^ < 0, > 0, they 
act as gaps for atoms and diatoms respectively. There is no Fermi surface in vacuum, hence fi^ < 0. The system is 
above criticality in the vacuum limit, i.e. it is in the symmetric phase; hence > 0. These two conditions define a 
quadrant in the — /x^ plane. Moreover, due to the non-relativistic nature of the problem, the zero energy level can 
be shifted by an arbitrary constant. This is a result of the symmetries of our models. The real-time {t = —it) version 
of the microscopic action S = Fa in coordinate space it, x) is symmetric with respect to the energy shift symmetry 

^__e*^V y'^e^'-^V n^fi + E. (7) 

Since no anomaly of this symmetry is expected and our cutoff respects this symmetry (see below), this is a symmetry 
of the flow equations and the effective action Tk=o- Hence, by the appropriate energy shift, we can make one energy 



^ Flavor indices applicable for Systems II and III are suppressed in this section. 
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state gapless, i.e. put it on the boundary of the quadrant in — /i^ plane. We end up with three distinct branches 

M- 

w? > 0, fi^ — atom phase (a^^ < 0), 

m| = 0, < diatom phase (a^^ > 0), (8) 

= 0, = resonance (a~^ = 0). 

For System II the gapless state is the lowest energy state. In the atom phase (a^^ < 0) diatoms (p are gapped and 
the lowest excitation is an atom tp. In the diatom phase (a^^ > 0) the situation is reversed: (p is the lowest excitation 
above the vacuum and tp has a gap — /x^, which can be interpreted as a half of the binding energy oi ip, e = 2^^. 
At resonance (a^^ = 0) both, ip and tp, are gapless. For Systems I and III, and for small values of |a|^^, one finds 
a whole spectrum of trions, bound states of three atoms, which have a lower energy then atoms and diatoms. This 
effect has been first predicted and calculated by Efimov in a quantum mechanical computation 4] , and modifies the 
vacuum structure j2l| . In the trion phase both atoms and diatoms show a gap, i.e. the ground state has fi^ < 0, 
771^ > 0. However, for an investigation of the excited Efimov states we may as well use the vacuum fixing condition 
([8|. At resonance this corresponds to degenerate energy levels of the Efimov states and the atoms/diatoms, which 
becomes a good approximation for the high Efimov states which are close to the atom/diatom threshold [2lj . 

The vacuum limit, which we described above, leads to numerous mathematical simplifications. For example, all 
diagrams with loop lines pointing in the same direction vanish in the vacuum limit. This can be demonstrated using 
the residue theorem for the frequency loop integration. Indeed, the inverse propagators have non-negative gaps and 
all considered diagrams have poles in the same half plane of the complex loop frequency. Thus we can close the 
contour such that it does not enclose any poles and the frequency integral vanishes. The argument works also for 
the IPI vertices provided they have poles in the same half-plane as the propagators. This finding simplifies the RG 
analysis in vacuum considerably. For example, one can show that in vacuum the atom inverse propagator P^{Q) is not 
renormalized |3^ . The only one-loop diagram, which renormalizes P^, has inner lines pointing in the same direction, 
and therefore vanishes. It is sufficient to analyze only one-loop diagrams because the RG flow equation ([T]) has a 
general one-loop form (23j . Another very important simplification in vacuum comes from a special hierarchy, which 
is respected by the flow equations. We define the n-body sector as a set of 2n-point IPI vertices written in terms of 
elementary atoms (in this sense P^{Q) belongs to the two-body sector because ip ~ ipip is composed of two atoms). 
The vacuum hierarchy consists in the fact that the flow of the n-body sector is not influenced by any higher-body 
sectors. The flow equations for the n-body sector simply decouple from the fiow of the (n+l)-body sector (and higher). 
The observed hierarchy is a consequence of the diagrammatic simplification in vacuum. At finite density (n ^ 0) or 
temperature (T 7^ 0) the decoupling of the low n vertices from the high n vertices is not valid anymore. 



IV. TWO-BODY SECTOR: EXACT SOLUTION IN THE DIATOM PHASE FOR A POSITIVE 

SCATTERING LENGTH 



The two-body sector truncation is defined by: 

Ffc = rfe(2) + rfe(3) (9) 

in all three models ([3]), (01) and ([5]). As mentioned in Sect. Illll the RG flows belonging to the two-body sector decouple 
from higher-body sectors in vacuum. Due to the non-renormalization of the atom propagator, it is sufficient to solve 
the fiow equations only for the Yukawa coupling h and the diatom inverse propagator P^. 
It turns out that the Yukawa coupling is not renormalized in vacuum for all three models: 

dth = 0. (10) 

Due to the U(l) phase symmetry, there is no one-loop Feynman diagram in our truncations (I3|4|5p . which renormalizes 
the Yukawa coupling h. The only one-loop diagram, which could contribute to the flow of h, contains the four-atom 
vertex A^((3i, Q2, Qa)- The vertex A^((5i, (52, Q3) is not renormalized in vacuum (see Appendix iBjl and vanishes on 
all scales, provided its microscopic value is zero. The argument can be extended to a momentum-dependent Yukawa 
coupling h{Qi,Q2,Q3)- 

In order to solve the two-body sector, it remains to calculate the flow of the diatom inverse propagator P^{Q), 
which is schematically shown in FIG. ([T]) and can be written as follows: 

dtP^iQ) = -j^ J dt ^p^^^^ ^ R^{L)){P^{Q -L) + R^iQ - L))' ^^^^ 

L 



a — . — = a 



■o- 



FIG. 1: Schematic graphical representation of the flow for the inverse diatom propagator P^. Diatoms are denoted by dashed 
lines, atoms by solid lines. 

where p = +1 for bosons and p = —1 for fermions. It turns out that the flow of the inverse diatom propagator in 
System III is exactly the same as in the System II. In the last formula one has: 

P^^iQ) = (itoq + q^ - t^^) (12) 

and R^{Q) stands for the atom regulator. 

It is remarkable that using a special choice of the regulator, we can integrate the flow pd|) exactly. We follow p^[33| 
by choosing a regulator, which is frequency and momentum independent and has the form i?0 = = h? exp 2t. This 
regulator has the advantage that it is Galilean invariant and hence the Galilean symmetry of the microscopic action 
is preserved during the RG evolution. First we perform the frequency loop integration in Eq. pTjl with help of the 
residue theorem: 

(r.. ^ [ ^ ^ 

* 3 + pJ (27r)3 *iWq + 12 + (l-q)2-2/x^ + 2i?^ 

2/i2 r (Pi 1 



3+p 7 (27r)3 iwq + P + (1 - q)2 - 2Ai^ + 2i?^ 

where the second equality holds due to the non-renormalization of the Yukawa coupling h and atom inverse propagator 
P^, and thus dt — > dt- Using the specific values of the regulator R^{tjii) = and Rjp{t = 0) = we integrate out 
the flow equation from the UV scale i = to the IR scale tjjf = — oo and obtain: 

.iR(n^ uuv^^^ 2/.2 r dH f I 1 



P'''{Q)-Pr{Q) = 




12 + (1 - q)2 - 2^^ iuj^ + 12 + (1 - q)2 - 2^1.4, + 2A2 

P /2 \ 



2 



(A-^/^ + ^-M.). (14) 



47r(3+p)^ V 2 4 



The last identity assumes A >> |^^|, |q|, jc^ql. 

At this point we must fix the initial condition P^^{Q) at fc = A in order to obtain the physical inverse propagator 
P^^iQ) at = 0. This is done in [l^, [sil, [12] and we follow the same steps here. For broad resonances with h? ^ 00 
the inverse diatom propagator at the microscopic scale A is given by: 

p!1^ ^ v{B) + Sv, iy{B) ^ fiB{B - Bo). (15) 

Here i'(-B) is the detuning of the magnetic field B which measures the distance to the Feshbach resonance located at 
Bq. The magnetic moment of the diatom is denoted hy fiB- The counter term Sv depends on the ultraviolet cutoff 
A. Neg lecting a possible background scattering length atg, the scattering length a and the detuning i^(i?) are related 
by mi: 

(16) 



471(3 +p)v{B)' 



Thus the Yukawa coupling is proportional to the square root of the width of the Feshbach resonance. For narrow 
Feshbach resonances (h — s- 0) perturbation theory is applicable, while for broad Feshbach resonances {h — > 00), which 
are of main interest in our work, the problem becomes strongly coupled. Using Eq. (|15|) and (jl6p . we rewrite Eq. 
(El) as: 
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At this point the momentum independent counter term 5v can be identified: 



= ^tJ^A (18) 
47r(3+p) ^ ^ 



and we obtain our final result for the fc-dependent inverse diatom propagator P^p^kiQ) 



2 



The wave-function renormalization Z^p^k can now be defined: 



d{iLUq 



/l2 1 



,=0 M^+P)j^/k^-fi^' 
z 



(20) 



and the IR inverse diatom propagator P^piQ) reads: 

In vacuum and for positive scattering length (a > 0) the vacuum condition, Pip{Q = 0) = = 0, must be 
fulfilled. This leads to: 

a=-^. (22) 

For positive scattering lengths in vacuum — /i^ is a positive gap of the atom tp and can be interpreted as half of the 
binding energy of the diatom e^. Hence, the binding energy can be expressed as: 



2 1 

Ma2 



= - -1772- (23) 



The second equation is expressed in conventional units and is the well-known universal relation for the binding energy 
of the shallow diatom Q . It should be mentioned here that the two-body sector can also be solved exactly using a 
non-relativistic version of the Litim cutoff , which is optimized in the sense of [H, [s^ . The drawback of this cutoff 
is that it breaks Galilean symmetry and one has to put some Galilean non-invariant counter terms into P^^ {Q) to 
restore Galilean symmetry in the IR. 

It is important to stress the appearance of universality in the broad resonance limit {h^ — > oo) (3lj |: The IR physics 
becomes insensitive to the initial conditions in the UV. For example, one may consider possible momentum-dependent 
modifications of the microscopic inverse propagator P^^k=K, which result in deviations from an exactly pointlike 
form. Their effect on P^^k=o is suppressed by with respect to the quantum loop contribution and it therefore 
becomes irrelevant in the broad resonance limit. In vacuum, and for h — )■ oo, the only physically relevant scale is 
given by the scattering length a. 

Physics becomes completely universal if we perform the unitarity limit, h? ^ oo (broad resonance limit) and 
— *■ (resonance limit) [37| . In vacuum, all scales drop out in this limit. The atom and diatom inverse propagators 
take the following form: 



P^{Q)=iu^ + ci\ P^^Q)^Z\j'-^ + ^. (24) 

An alternative quantum-mechanical derivation of Eq. (|24p can be found in Appendix O 

Let us perform a scaling dimension counting in the unitary limit^. We start with the fact that [F] = 0. In 
non-relativistic physics, energy scales as two powers of momentum^ and the free field scaling reads: 

[H = 2, [V']=3/2, [^]=3/2, [h] = l/2. (25) 



^ We denote a scaling dimension of some quantity X by [X]. 
* This is known as the dynamical exponent 2 = 2. 
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For the universal interacting theory the scaHng of is modified according to Eq. p4)) . The scahng dimension of the 
diatom field and Yukawa coupling h can be computed from the Yukawa term and the kinetic term of (p: 

\h\ + \ip\ + 2[V^] = 5 Yukawa, 2[/i] + 2[(^] + 1 = 5 Kinetic. (26) 

This system is degenerate and we get a solution [/i] = a and [(^] = 2 — a, where a is some real number. The absence 
of a scaling of h in Eq. ((TO)) . however, fixes \h\ = and [yj] =2. Note that the scaling of the diatom field at unitarity 
is different to the scaling of the atom field i/'- This is a manifestation of the fact that the scaling in the two-body 
sector is governed by a fixed point'^ , which is different from the Gaussian fixed point [s^, [sH] . Exactly at unitarity no 
obvious scales are left in the problem and the theory seems to be scale invariant. Even more, at the two-body sector 
level the theory seems to be an example of a non-relativistic conformal field theory (NRCFT)^. This type of theories 
are symmetric with respect to the Schrddinger group, which is an extension of the Galilean symmetry group (for more 
details see Appendix |A). It is known, however, that the Schrodinger symmetry can be broken by a quantum anomaly 
in higher-body sectors 5|. The fate of the Schrodinger symmetry is different for the different systems considered. For 
the resonantly interacting particles (Systems I and III) it was demonstrated by Efimov [3| a long time ago that in 
the three-body sector the continuous scaling symmetry, which is a part of the Schrodinger symmetry, is broken to the 
discrete scaling subgroup Z |3j] . This manifests itself in the appearance of a geometric spectrum of bound states in 
the three body sector, which is called the Efimov effect. For the System II of SU(2) symmetric fermions it is believed 
that the Schrodinger symmetry is not broken in the higher sectors of the theory and that this is a real example of an 
NRCFT [SQ^. 

To summarize, in this section we have solved exactly the two-body sector in vacuum for a positive scattering length. 
The solution (|2ip was obtained for the specific initial conditions hk=A{Qi,Q2) = h, X^^k=A = and Pip^k=K{Q) 
given by Eq. (jlSp . This choice corresponds to a pointlike microscopic atom interaction. However, the presented 
calculations can be generalized to an arbitrary boson mediated atom interaction with A^,fc=A = while /ife=A(Qi, Q2) 
and Pip,k=K[Q) can be chosen freely. 



V. THREE-BODY SECTOR: FLOW EQUATIONS 

The main emphasis of this work is devoted to the analysis of the three-body sector of the three models (U) 
and ([5]) in the unitarity limit. We demonstrate the appearance of the Efimov effect in Systems I and III and its 
absence in System II from the field theoretical RG perspective. In the present section we formulate a flow equation 
for the coupling \d.{QijQ2,Q3) and make some general simplifications. In the next section we use the pointlike 
approximation for A3((3i, Q2, Qa) to solve the problem. The last two sections are devoted to the solution of the 
general momentum-dependent form of the flow equation. 

The closed, exact solution for the two-body sector provides a simple strategy for a computation of the coupling A3 
in the three-body sector. In general, one may introduce separate cutoffs and R^p for the atoms "0 and diatoms 
ip. The presence of the cutoff R^p does not affect our computation in the two-body sector. We may therefore first 
lower the cutoff R^ from to zero, while keeping R^ fixed, and subsequently lower R^ to zero in a second step [23j . 
As the result of the first step the diatom inverse propagator is modified according to Eq. (PT|) . This step also 
induces diatom interactions, as for example a term ~ ((^*(/3)^. However, these interactions belong to the four-body 
and higher sectors. By virtue of the vacuum hierachy, they do not influence the flow of A3. For the second step of our 
computation we can therefore use a version of the flow equation where only the diatom cutoff R^p is present. In this 
flow equation and P^p are flxed according to Eqs. (pij) and 

For the diatoms we use a sharp cutoff: 

i?.(Q,/c) = P,(Q)(^^-l). (27) 
The special feature of this cutoff is that the regularized diatom propagator takes a simple form: 

P^{Q) + R^{Q,k)^^^^'^^-^'^PM- ^^^^ 



^ called a unitarity fixed point 

^ Another example of NRCFT in two spatial dimensions is a theory of anyons [39l |. 
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FIG. 2: Kinematics of the vertex X3{Qf,Q2,E) in the center-of-mass frame. The atoms and diatoms have momenta Q 
(-E^i, qi) , Qi = {~E^i + E, -qi) and = (£^2, q2), Q2 = {-E^2 + E, -q2). 



Model 


A 


B 


C 


System I 


1 


2 


1 


System II 


1 


-2 


1 


System Ill(a) 


1 


-2 


1 


System Ill(b) 


4 


4 


1 



TABLE I: Numerical coefficients A, B and C in the flow equation (|29|l for the three examined systems. In the case of System 
Ill(a) we consider the scattering of the type V'iVi ~* ^ifj with i 7^ j, while System Ill(b) corresponds to the the vertex 

A3 = 3A3a + A36 



Thus the propagator is cut off sharply at the sliding scale k. Our choice of the cutoff is motivated by technical 
simplicity as well as effective theory [sj and quantum mechanical Q approaches to this problem. The advantage of 
this cutoff is the property of locality in the momentum space, which means that it chops off momentum shells locally. 
In the three-body sector we are interested not only in the IR value of the atom-diatom vertex A3 , but also in the flow 
at all scales. 

Let us now calculate the flow equation of the fPI atom-diatom vertex A3. For SU(3) fermions there are two atom- 
diatom vertices, A3a and As;,, and we postpone the analysis of this model to the end of the section. In Minkowski 
space (the real time version of our theory) the atom-diatom scattering amplitude is given by the amputated connected 
part of the Green's function (O|iy9?/'(y3^'0^|O), and thus it can be simply calculated from the knowledge of A3. We first 
consider the kinematics of the problem. The IPI atom-diatom vertex \^{Qi^ Q21 Qa) depends generally on three four 
vectors, i.e. six independent rotation invariant variables in the center-of-mass frame. We take the incoming atom 
and diatom to have momenta qi and — qi, and energies E^i and E — E^i, while the outgoing atom and diatom 
have momenta q2 and — q2 and energies E^2 and E — E^2- We denote the vertex in the center-of-mass frame by 
hiQt.Qt.E) (see FIG 2.). This configuration is in general off-shell which is necessary since, in the flow equations, 
the vertex also appears inside a loop. 

In Minkowski space^ the flow equation for the atom-diatom vertex A3 for the Systems I and II reads: 



dtX3{Qt,Qt;E) 



dt 



\l\-k) 



P^{L)P^{-L + Q) 



CX3iQt,L;E)X3{L,Qt;E) 



(29) 



A 



X3iL,Qt;E) + X3iQt,L;E) 



P4,{-L + Ql)' 



P4-L + Qf) P^{-L + Ql) 

where Q^Q^ + Qf = {E, 0). The coefficients A, B and C for Systems I, II can be found in TABLE I. The graphical 
representation of this equation is depicted in FIG. [3l The scale derivative on the RHS acts only on the cutoff and can 
be computed easily, (3t6'(|l| — k) = —kd{\l\ — k). 



The flow equation of the effective action JTJ is formulated in Euclidean spacetimc (imaginary time formalism). In order to obtain 
the flow equation in Minkowski space, it is sufScient to take external frequencies iOext to be imaginary, i.e. perform a transformation 
^a-^t — > ii^cxt, which is the inverse Wick rotation. 
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Fortunately, the flow equation can be simplified considerably. First note that there is only one inverse propagator 
P^{L) with a loop momentum L of positive sign in Eq. (P^ . For this reason the whole integrand in Eq. has 
a single frequency pole in the upper half plane. Thus the frequency integration in Eq. ([29|) can be performed with 
the help of the residue theorem by performing the substitution uj\ — > . This puts the atom in the loop on-shell, 
corresponding to Pti,{L) in Eq. (|29p . Wc obtain a simpler equation if we also put the energies of the incoming and 
outgoing atoms on-shell {Qf = (iqi,qi), ~ (^q^jQa))- The diatoms in the loop in Eq. ([^5]) are generally off-shell. 
To solve this "half-off-shell" equation only the values A3(qi, q2, E) = X^iQ^ — (iq^, qi), — (iq^, q2): E) are needed 
i- 

Our aim is the calculation of the atom-diatom scattering amplitude at low energies and momenta. For low momenta 
the dominant contribution is given by s-wave scattering. In principle, the right hand side of Eq. (|29p has also 
contributions from higher partial waves, which we neglect in our approximation and simplify the flow equation (|29p 
by projecting on the s-wave. This is done by averaging Eq. (|29p over the cosine of the angle between incoming 
momentum qi and outgoing momentum q2. Introducing the averaged IPI renormalized vertex, which depends on 
three scalar variables: 



1 

X3{qi,q2,E) = — J d{cos9)X3{qi,(i2,E), 



(30) 



we end up with the flow equation: 



dt\3iqi,q2,E) = - ^^^ + [CAgfa, fc, £;)A3(fc, g2, 



4 2 



I {A3(gi, fc, E)G{k, + G{qi, fc)A3(fc, 52, E)} + AG(gi, fc)G(fc, q^) 



(31) 



where the symmetric function G{qi, 92) is defined by®: 



G(9i, 92) = log 2 , 2 1 (32) 

49192 9? + 95 - 9192 - t 



The infinitesimally positive ie term arises from the Wick rotation and makes both Eq. (|3T|) and ([32)1 well-defined. It 
is remarkable that Eq. (|3ip is completely independent of the Yukawa coupling h and thus is a well-defined equation 
in the limit of infinite h. 

For SU (3) fcrmions the situation is more complicated because there are two vertices X^a and X^b h^ our truncation 
To extract the flow equation for A3Q we consider the scattering channel ipitpi — > ^j'fpj with i ^ j (e.g. ipitpi — > 
f2'4'2)- After performing the same steps as for System I and II, we end up with a flow equation: 

2('3 -I- ri) fc3 

5tA3a(9l,92,^^) = ^ , =[3A3a(9l,fc,£:)A3a(fc,g2,£^)+2A3a(9l,fc,g)A3b(fc,g2,£^) + 

2 {A3a(9i, k, E)G{k, (72) + G(gi, k)X3a{k, 92, £^)} + 

+ {A36(9l, fc, E)G{k, (72) + G{qi, fc)A3b(fc, q2,E)} + G{quk)G{k, 52)] , (33) 

where p = —1 and G{qi,q2) is defined in Eq. ([5^ . Note, that the coupling A3b appears in the flow equation for 
A3a. The flow equation for A3b can be extracted by considering the scattering channel (piipj ^iipj with i j (e.g. 

dtX3biqi,q2,E) = - ^^^+P^ ''^ [X3biqi,k,E)X3b{k,q2,E)+ 

- {A3b(9i, A:, E)G{k, 92) + G{qi,k)X3b{k, q2, E)} + G{qi, k)G{k, 92)] . (34) 

This equation is completely decoupled from Eq. ([55]) and has exactly the same form as Eq. (|3ip for SU{2) fermions. 
The reason for this is simple: The RG equation (l34|) has the graphical representation depicted in FIG. [S] It turns out 



* This is in fact a s-wave projected tree (one-particle-reducible) contribution to the fully connected atom-diatom vertex X3{qi,q2, E). 
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FIG. 3: Graphical representation of the flow equation for A3. Full lines denote atoms ip and dashed lines diatoms ip. The 
shaded circle denotes A3. 



that in this channel only one type of diatom (in our example ^2) and two types of atoms (■(/'i and ip^) appear, which 
is exactly the same as in the case of SU{2) fermions. Remarkably, it is possible to introduce a linear combination 
A3 = SAsa + Asf, for S'[/(3) fermions, which has a simple flow equation of the form (PT|) with coefficients A, B and C 
given in TABLE I (forth line). We call this System Ill(b). In the SU(3) fermion model the diatom-atom in-state 
(fiiipi can lead to the different diatom-atom out-states (ySiV'ij ^2^2 and (fz^'a- If the diatom-atom out-state is not 
a final but only an intermediate state(e.g. one is interested in the scattering into a three atom final state), we 
must sum the scattering amplitudes for all possible atom-diatom pairs. It easy to show that A3 = 3A3a + A3fc cor- 
responds to the IPI contribution to the full scattering amplitude ipiipi — > anything (e.g. ipiipi <y5iV-'i + '/'2'02 + <y53''/'3)- 

To summarize, although at first sight it seems that for SU(3) fermions we must solve a system of two flow equations, 
it turns out that for the two specific situations it is sufficient to solve only one equation pip . This equation is the 
main result of this section. In the next sections we solve this final version of the RG flow equation for atom-diatom 
IPI vertex for all three systems using various approaches. 



VI. THREE-BODY SECTOR: POINTLIKE APPROXIMATION 



In this section the flow equation ([5T|) will be solved employing a simple and intuitive pointlike approximation. The 
IPI vertex X3{qi,q2T E) will be replaced by a single momentum- independent coupling A3(£'). In the low energy limit 
(E — !■ 0) the flow equation pip takes a simple form in the pointlike approximation: 



j + fA? + C(Af)^' 



2Af, (35) 



where we use G{q —^Q,k)-^ from Eq. ([5^ . The renormalized coupling is defined as A|^ — X^k^. This definition 



2k^ 



is motivated by a simple power counting near the unitary fixed point {[iptpip*ip*] = 7 — > [A3] = —2 — > [A3 ] — 0). The 
RHS of Eq. psp is a quadratic polynomial in A|* with constant coefficients. This type of equation is discussed in 
Appendix [D] The behavior of the solution is governed by the sign of the discriminant D of the quadratic polynomial 
on the RHS of Eq. ([35]): 



• D > ~ fixed point solution 

• = 0- see Appendix ID] 



D < ~ periodic limit cycle solutions with a period T = 



2ir 



-D ' 



The discriminant is given by: 



In the special case of the Systems I, II and III the solution in the pointlike approximation is summarized in TABLE 
II. For Systems II and Ill(a) we find the solution with a fixed point with vanishing anomalous dimension 77 = in 
the IR (see Appendix [Pjl . For Systems I and Ill(b) the situation is completely different. Wc obtain a periodic limit 
cycle solution of the form X§{t) tan [Tt]. The intuitive interpretation of this solution is that during the RG flow 
we hit three-body diatom-atom bound states, which manifest themselves as divergences of A^. In the unitary limit 
there are infinitely many of these bound states, which are equidistant in a logarithmic scale. The continuous scaling 
symmetry is broken to the discrete Z group. This is the well-known Efimov effect 0,13 1 which indeed is present for 
equivalent bosons (System I) and is absent in the case of SU(2) fermions. In the case of equivalent bosons (System I) 
the Efimov result is: 

^ = exp(-27r/so) (37) 
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Model 


D 


r 


so 


System I 


-7.762 


2.255 


1.393 


System II 


9.881 






System Ill(a) 


9.881 






System Ill(b) 


-7.762 


2.255 


1.393 



TABLE II: Discriminant D, temporal RG period T (if applicable) and Efimov parameter so (if applicable) in the pointlike 
approximation for Systems I, II, Ill(a) and Ill(b). 



with En+i and £"„ denoting neighboring bound state energies. The Efimov parameter sq is given by the solution of a 
transcendental equation and one finds sq ~ l.OOGM. By dimensional arguments we can connect the artificial sliding 
scale fc^ with the scattering energy E as E ^ k"^ [21[ . The proportionality factor disappears in the ratio of the energies 
and hence the Efimov parameter can be read off from the RG period: 



fc2 - - exp(-2r) ^ so - -. (38) 

The values of the Efimov parameter for Systems I and ni(b) can be found in TABLE II. We obtain so ~ 1.393, 
which differs from the correct result by 40%. In the next two sections we demonstrate that the simple pointlike 
approximation is too crude to get the correct quantitative agreement. Nevertheless it provides us with the first hint 
how the Efimov effect appears also in the functional renormalization group framework. 



VII. THREE-BODY SECTOR: SYSTEMS I AND II 



In this section we only discuss Systems I and II leaving the analysis of System III to the next section. It turns out 
that in these two cases the flow equation ((3T1) for ii^ = can be formally solved exactly. For two-component ferniions 
this was shown by Diehl et. al. in [Tsj . To find the exact solution most easily we perform the following redefinition: 

ft{hM,E) = A{i+p)qiq2\z{qi,q2.E) g(ti, ^2) = 4(3 + p)gi<Z2G(gi, 92), (39) 

where, from now on, we prefer to work with logarithms of momenta ti — ln(qi/A) and ^2 — In(q2/A). As before 
p = -fl(p = — l)in the case of bosons (fermions). The RG scale dependence of the reduced atom-diatom vertex 
ft{ti,t2, E) is denoted by the subscript t. It is important to stress that we are generally interested in the solution 
of Eq. ([3T|) for the scattering of particles of non-zero energy E. Nevertheless, we observe that the energy E cuts off 
the RG flow in Eq. (|3ip in a similar way as the regulator ([?f|) . With this relation between and E in mind, the 
coupling for fc 7^ and E ^ imitates the effect of a non-zero energy of the scattering particles, i.e. k ^ 0, E ^ Q. 
The flow equation at vanishing energy E = now reads: 



dtftitiM) ^ 



Ag{ti,t)g{t, + 1 {ft{ti,t)g{t, + 9{ti,t)ft{t, 12)} + C ft{h,t)ft{t, 12) 



(40) 



We assume that in the UV the reduced atom-diatom IPI vertex is vanishing, i.e. the initial condition is /t=o(^i; ^2) — 0- 
In general we are dealing with the Riccati differential equation in matrix form, where both matrices g and ft have 
a continuous index running in the interval ^1,^2 G (— oo,0). The RHS of Eq. is a complete square, which is a 
special feature of the Systems I and II (see TABLE I). In order to flnd the formal solution of Eq. (|40|) we define: 

Utut2)^pMtut2)+g{h,t2), (41) 

which can be recognized as the reduced, fully connected atom-diatom vertex. The fiow equation for the full vertex 
/t(^ii^2) with the initial condition takes the simple form: 

dtft{tiM) = -^ft{ti,t)!t{t,t2). ft=o{tiM) = 9{tiM)- (42) 

It is convenient to rewrite Eq. (|42p in matrix notation (/t(ti,t2) ~^ ft)'- 

dtft = -^ft ■ At ■ ft, ft=o = 9, (43) 
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where At has matrix elements At(ti, 12) = 5(t—ti)d{t—t2) and matrix multipHcation denotes t-integration. Multiplying 
both sides of Eq. from the left and right by ff^ we obtain: 



dtft"' = -ft~'-dtffft^' = ^At, ft-Jo=9~\ (44) 

VoTT 



which is formally solved by: 



ft^ [1+^ I dsg-A,] -g (45) 





for t G (— cx),0). / denotes the identity matrix. 

In the IR limit t — > —00, which corresponds to integration of all quantum fluctuations, / = ft=-oo solves the 
following matrix equation: 

J^g+^g.f. (46) 

This is the well-known STM integral equation for bosons {p — +1) and fermions (p = —1) for the half-off-shell, 
amputated, connected Greens function^ Q. 

The difference of the signs in Eq. ([IS]) between System I and II turns out to be crucial. In order to see that, we 
solve Eq. ([iS]) numerically by discretization. A series of cartoons of the evolution of the reduced IPI vertex ft{ti,t2) 
for both systems is shown in FIG. SI For fermions, first a peak appears in the UV {ti = 0,t2 = 0), which propagates 
in the diagonal direction {ti — ^2) during the RG evolution. On the other hand, for bosons, a periodic structure (with 
period Tspatiai ~ 6.2 in both directions) develops gradually. Now it is clear why the approximation investigated in 
the last section failed to give the quantitatively correct result. The pointlike approximation, which corresponds to a 
planar landscape (no ti and t2 dependence, see Section IVl)) . is not valid in the three-body sector (for more details 
see Appendix [U]) . The evolution in the RG time t of the UV point /t(ii = 0,t2 = 0) for both systems is depicted 
in FIG. [5l While for fermions the evolution is monotonic in time, in the case of bosons we obtain a "temporal" 
oscillation of period Ttemp ~ 3.1. For different points in the ti — 12 plane the "time" evolution is triggered at the scale 

tin ~ 0{ti,t2)- 

The numerical solution for bosons is consistent with the results of [3, Q. Spatial and temporal oscillations are 
correlated. As found in [1, Q evolution in the RG time develops a limit cycle behavior. The Efimov parameter sq can 
be calculated Sq = — « 1.0 (see Section IVH) which is in a good agreement with the Efimov result sq « 1.00624. 
The accuracy of our result is limited by the numerical procedure only. 



VIII. THREE-BODY SECTOR: SYSTEM III 



As introduced in Sect. |V]for SU{3) fermions there are two specific situations (System Ill(a) and System 111(b)), 
when there is a single flow equation instead of the general two. Fortunately, both cases can formally be solved for 

= in a similar fashion compared to Sect. IVIII In fact. System Ill(a) is completely equivalent to System II (see 
TABLE I) such that we obtain a flxed point solution in this case (see FIG. [5]). For System Ill(a) we follow similar 
steps as in Sect. IVIII we define a reduced atom-diatom IPI vertex ft{ti,t2, E) ((39| and obtain a flow equation for 
the reduced vertex with the coefficients A — A, B = 4 and C ~ 1. These coefficients form a complete square and 
thus it is useful to define the fully connected atom-diatom vertex: 

UhM) ^ MhM) + MtiM)- (47) 

The flow equation ([^0]) now reads: 

dJt{h,t2)^-^Uti,t)ft{tM), ft=o(.ti,t2)=2g{h,t2). (48) 

The equation and the initial condition are identical to Eq. (|^^ for bosons^". For this reason we expect the appearance 
of the Efimov effect for the SU{3) fcrmionic System Ill(b) with the Efimov parameter sq ^ 1.00624. 



^ up to our redefinition II39I I 

^ The initial condition for SU (3) fermions is /t=o(<l , <2) = 2g(ti , t2), while for bosons one has ft—oitl tt2) = 9{tl,t2)- However, for bosons 
ff(ilii2) is two times larger than for fermions II39I I. Thus the initial conditions are identical. 
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FIG. 4: The RG evolution of the momentum dependent modified vertex /t(ti, t2) ~ 4:{3 + p)qiq2X3{qi, q2, E) for bosons (al-a4) 
and SU(2) fermions (bl-b4). Spatial momenta ti, t2 and the RG time t are descritized to A'^ = 50 intervals with a step At = 0.4. 
The cartoons for bosons and fermions correspond to the descritized steps 10, 25, 35, 50. 

At first sight it seems surprising that both bosons and SU{2>) fermions have the identical Efimov parameter sq. As 
explanation, we propose a simple possible quantum mechanical argument: In order to find a bound state spectrum 
for SU(3) fermions one must solve the three-body Schrodinger equation. The total wave function must be totally 
antisymmetric for fermions. We can achieve this by taking the total wave function as the product of a totally 
antisymmetric flavor part (eyfe|i > \j > \k >) times a totally symmetric orbital part. Hence the orbital part has the 
same symmetry property as the bosonic case. Only the orbital part is needed for the quantum mechanical calculation 
of the bound state problem, which leads to the identical Efimov parameters for bosons and SU(3) fermions. 
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FIG. 5: Numerical evolution in the RG time t of /t(ti = 0,f2 = 0) for System I (A) and System II (B). For SU(2) fermions 
(B) the modified vertex approaches a fixed point solution; in the case of bosons (A), a limit cycle behavior is developed with a 
period Ttemp ~ 3.1. 



IX. CONCLUSIONS AND OUTLOOK 



This work applies the method of functional renormalization to the few-body physics of atoms near a Feshbach 
resonance. We investigate three different systems, namely identical bosons as well as two and three species of fermions. 
The two-body sector is solved exactly. The unitarity limit is governed by a fixed point and all three systems seem 
to be examples of the non-relativistic conformal field theories. In the three-body sector, however, no infrared fixed 
point exists for bosons and three-component fermions. We solve the momentum-dependent problem of the three-body 
sector at unitarity. This leads to the Skorniakov-Ter-Martirosian equation, well-known from quantum mechanics. A 
numerical solution for U{1) bosons and SU{i) fermions shows the emergence of the Efimov effect; the appearance of 
an infinite geometric spectrum of triatom states. Hence in these systems the continuous scaling symmetry is broken 
to the discrete scaling subgroup Z by a quantum anomaly. The renormalization group flow develops a limit cycle 
behavior (see FIG. [5]). The Efimov parameter sq for the three-component fermions is found to be identical to the 
Efimov parameter of the well-studied bosonic case, which agrees with the quantum-mechanical prediction. 

The current work can be extended in various ways: One can go away from unitarity in the three-body sector and 
derive universal properties such as recombination rates and the positions of diatom-triatom thresholds. Our technique 
allows us to investigate equilibrium states with non-zero density and temperature. This can be achieved by simply 
changing the chemical potential and introducing the temperature by replacing the w-integrals by the discrete sums 
of the Matsubara formalism. In that case the effective, approximate description of the models in terms of various 
composite fields (e.g. trions, density bosons) would be very useful due to a large reduction of the numerical effort. 
The description of a simple, but efficient effective theory is summarized in Appendix [E] The excellent agreement of 
the vacuum solution with high precision quantum-mechanical computations provides a robust starting point for the 
investigation of the many-body system of nonzero density and temperature. 
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APPENDIX A: GALILEAN AND NON-RELATIVISTIC CONFORMAL SYMMETRY 



All systems we consider in this work ([3|), (|4]) and ([5]) have a centrally extended Galilean spacetime symmetry^^. 
The centrally extended Galilean algebra consists of eleven generators: particle number TV (central charge), time 
translation 77, three spatial translations Pi, three spatial rotations My and three Galilean boosts Ki. The non-trivial 
commutators are (in the real time formalism): 

[My , Mki] = i{5ikMji - SjkMu + SuMkj - SjiMkt), (Al) 



This is the non-relativistic analogue of the Poincare group in relativistic QFT. 
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[M,„ i^fe] = HS^kKj - 5,kKi), [M,,, Pfe] = i{5,kP, - 5,kPi), (A2) 

[P„ K,] = ^id^jN, [H, K,] = -zPj. (A3) 

In the case of a free non-relativistic field theory the group of spacetime symmetries is in fact larger than the Galilean 
group [i^ |4l| and is called the Schrodinger group^^. For the dynamical exponent z = 2 there are two additional 
generators: the scaling generator D and the special conformal generator C. The scale symmetry acts on the time and 
spatial coordinates according to: 

{x„t)^{xlt')^{Xx,,Xh), (A4) 
where A is a scale parameter. A special conformal transformation on time and spatial coordinates is given by (40| : 

(x„t)-(x'„t')-(Y^,Y:^), (A5) 

where c is a parameter of the special conformal transformation. The additional, non-trivial commutators of the 
Schrodinger algebra are: 

[P„D] = -iP,, [P,,C]^^iK,, [K„D]^iK,, (A6) 

[D, C] = -2iC, [D, H] = 2iH, [C, H] = iD. (A7) 

It is important to note that besides the free theory there are few known examples of interacting theories which are 
symmetric with respect to the Schrodinger group. These theories are called non-relativistic conformal field theories 
(NRCFT) and SU (2) non-relativistic fermions at unitarity (System II) are believed to constitute one of them. 

In analogy to relativistic conformal field theories it is possible to introduce primary operators in an NRCFT [soj . 
A local primary operator ff(t,x.) has a well defined scaling dimension Aff and particle number Ng-: 

[D, ff] = iAffff, [N, e\ = Ne6, (A8) 

where ^ =#"(;; = 0, x = 0). The primary operator 6 also commutes with Ki and C: 

[K,, ^] = [C, e\ = 0. (A9) 

It is possible to show that the operators, constructed by taking spatial and time derivatives of a primary operator 
form an irreducible representation of the Schrodinger group. Similar to the relativistic case the form of the two- body 
Greens function of the primary operators is fixed by the conformal symmetry (in the imaginary time formalism) [39| : 

<^^t,.(,^_,_S_)^, (AlO) 

where v — A^y — 5/2 for d = 3. The simplest examples of primary operators in the theory of SU(2) symmetric fermions 
are the atom operator ?/' — 1, = 3/2) and the diatom operator ip {N^ — 2, — 2). The form of the inverse 
propagators at unitarity, which we found to be given by Eq. ((^ . is consistent with Eq. (jAlOp . 

APPENDIX B: COMPLETION OF THE VERTEX EXPANSION TO v'^^ 

In this appendix we complete the vertex expansion to fourth order and argue that our truncations pi4l5p are 
sufhcient to perform exact calculations for the three-body sector. At fourth order in the fields there are only two more 
vertices, which are compatible with the internal symmetries of the considered models: 

r(4)^ = i / \i.{Qi,Q2,QiWi.QM{Q2)^KQi)^{Qi)^{-Qi + Q2-Q-i + Qi). 

Qi,...Q4 

Ql:...Q4 



This is the non-relativistic counterpart of the conformal group. 
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In the two-channel model considered in this work, we choose the initial UV value of the vertex to be zero, = 0. 
This means that the interaction between the atoms are described at the microscopic level by the exchange of diatom 
states. With A^, = at the UV scale this coupling is not regenerated by the flow in vacuum. The one-loop diagrams 
contributing to the flow have inner lines pointing in the same direction with respect to the loop momentum and 
therefore vanish in vacuum (sec Sect. Thus A^ = is a fixed point. The flow away from this fixed point has 

been studied for the system with two species of fermions in [s^ ]. 

The IPI vertex A;^ belongs to the four-body sector (for definition of the n-body sector see Sect. IIIII) and it decouples 
from the flow equations of the two and three-body sectors due to the vacuum hierarchy (for more details see Sect. 
IIIII) . Thus our truncations (|3l4l5p are sufficient to obtain the exact vacuum physics of the three-body sector. 



APPENDIX C: BOUND STATE APPROXIMATION AND SEPARABLE POTENTIAL 

In this appendix we present an alternative solution of the two-body sector using the Lippmann-Schwinger equation 
of quantum mechanics, which helps to elucidate the efficiency of the two-channel model and the limitations of the 
trion approximation in [211] . 

The one-channel model provides an alternative description of ultracold atoms near a broad Feshbach resonance. 
This model contains the atom field 'ip only and the microscopic action is given by^^: 

rt=o = J r{Q){^uj + q^)m) + Y J i^*iQiMQ2WiQ3mQ4)6{-Qi + Q2-Q3 + Q4), (ci) 

where A^ is a pointlike four-atom interaction which is related to the s-wave scattering length in the IR. Roughly 
speaking, the quantum-mechanical atom-atom interaction potential of the one-channel model (jCip in Minkowski 
space is given by^**: 

y(x) = ^^(3)(x). (C2) 
Let us now perform a Fourier transformation of this potential: 

V{k,k')^ J dVexphi(k'-k) •r]T/(r) = ^. (C3) 

At this point two important remarks about the potential (jC3|) can be made: 

• V^(k, k') is a separable potential because it can be written in the form ^t/(k)t/(k'). 

• F(k,k') is k and k' independent, i.e. U{k) = U{k') = 1. 

We investigate the atom-atom scattering in the centcr-of-mass frame. The Lippmann-Schwinger integral equation for 
the K-matrix is [4^ : 

m.',^)^V„<,.-) + ^/^^2]i|^^, ,C4) 

where denotes the Cauchy principle value and A is a momentum cutoff, which regularizes the contact interaction. 
The K-matrix is similar to the T-matrix but uses a standing wave boundary condition which leads to the principal 
value prescription in Eq. (|C4p . The kinematics of K{k,k' , E) is similar to the kinematics depicted in FIG. [D The 
integral equation (|C4p can be easily solved in the special case of a separable potential. The solution factorizes: 

K(UX,E)^-mip, (C5) 



For simplicity, we present the one-channel model for U{1) bosons at unitarity only. However, our arguments can be extended away from 
unitarity and are applicable to both 5(7(2) and SU{3) fermion systems. 

Strictly speaking, the contact interaction is ill-defined and must be regularized. This can be done by introducing the pseudo-potential 
V{r)ip{r) = X^S^^'^ {r)-^{rip{r)). Here we use an alternative regularization by introducing a momentum cutoff A directly into the 
Lippmann-Schwinger equation. 



18 



where D{E) is given by: 



In the special case of the contact interaction Eq. (jCSp depends only on E. This means that the exact atom-atom 
scattering amplitude in the center-of-mass frame can be rewritten in terms of the exchange of the composite diatom 
with inverse propagator P,j){E, p = 0) ~ D{E). For the contact interaction D{E) is given by: 

^ dq <f 2 A E _^r^ dq 



The microscopic can be adjusted such that 



2 A 

A^. 47r2 



(C8) 



At unitarity the first two terms in the second Eq. (|C7|) cancel. The last integral is convergent, hence we take A ^ oo. 
By dimensional analysis we obtain: 



D{E) - VE. (C9) 

To summarize, the atom-atom scattering amplitude is momentum-independent in the center-of-mass frame. Hence 
the two-body sector can be solved exactly by introducing a diatom exchange in the s-channel. By Galilean symmetry 
this result can be extended to a general reference frame: 



The functional form of the inverse diatom propagator is consistent with our findings (|24p in Sect. IIVI 

In the three-body channel the atom-diatom interacting potential is momentum dependent in the center-of-mass 
frame. The momentum dependence is generated by the box diagram (see first diagram on RHS of FIG. [3]). For this 
reason the trion approximation, which we used in [2l| . and in particular the pointlike approximation (Sect. IVip . do 
not fully capture this momentum dependence and is not as efficient as the "diatom trick" . It leads to the quantitative 
inaccuracy of the Efimov parameter so- 



APPENDIX D: ANALYSIS OF ^/(t) = af{tf + /3f{t) +j 

In this appendix we perform an analysis of the differential equation, which we encountered in the calculation of the 
three-body sector in the pointlike approximation: 

^/(t) = a/(t)2 + /3/(t) + 7, /(to) = /, (Dl) 

with a, 7 e M. The form of the solution is determined by the sign of the discriminant of the /3-function D = — 
There are three different cases (without loss of generality we take a < 0, which is the case for Systems I, II and III): 

• £1 > 

In this case the /3-function has two fixed points /i (IR stable) and /2 (IR unstable) with /i < /2 (sec FIG. [HI). 
For the initial condition / < /i the solution is attracted to the fixed point /i in the IR. If the initial condition 
is / > /2, the solution is repelled from the fixed point /2 and in driven towards a Landau pole. Formally, we 
can extend the solution beyond the Landau pole. For an initial condition /i < / < /2 the solution is attracted 
by the IR fixed point /i and has the form: 



-13 + Vl3tanh \^{t + k) 
2^ 

with K fixed by the initial condition 



m ^ ^, (D2) 
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• L> = 

In this case there is one fixed point and we have four distinct subcases: 

a) q: = /3 = 7 = 0=>A trivial solution f{t) = / is obtained. 

b) a = /3 = and 7 7^ ^ We have a linear solution /(<) = / + j{t — to). 

c) a 7^ and /3 = 7 = ^ We obtain the Landau pole solution: 

fit) - .y (D3) 

1 + fa{t - to) 

This solution can be formally extended beyond the Landau pole. 

d) a ^ 0, (3 ^ and 7 7^ ^ This case can be put in the form of the case c) by a transformation / /+ •v/J- 

• D <0 

In this case there are no fixed points. The formal solution can be written as: 



-/3 + V^tan ^^^(t + n) 

fit) = (D4) 

Za 



where k is fixed by the initial condition. This solution is periodic with a period T = 



2tt 



-D ■ 



It is important to mention that the discriminant D is invariant under the multiplicative reparametrization f{t) — s- 
Lof{t) with some constant uj. This means that both the classification and the period T are not sensitive to the 
multiplicative reparametrization of in Sect. IVII 



APPENDIX E: BOSONIZATION AND FERMIONIZATION 



It was shown in Sects. I VIII and IVIIII that for a quantitative precise description of the three-body physics the 
momentum dependence of the vertex A3 is quite important. Although it was eventually possible to take this momentum 
dependence into account, we had to pay a price for it. Solving the flow equation numerically while taking the full 
s-wave projected momentum dependence of the vertices into account needs a relatively large numerical effort. While 
this is still manageable in the vacuum where both the density and the temperature vanish, the numerical cost would 
be significantly larger in the more general case of nonzero density or temperature. It is therefore reasonable to look 
for an effective approximate description that is numerically less expensive but nevertheless leads to good numerical 
precision. 

For this task it is crucial to find a simple way to take at least the qualitative features of the momentum dependence 
into account. For simplicity we concentrate this discussion on System III only. Let us start with the discussion of Asfc. 
From Eq. ^ we can read off that it corresponds to a channel where the fermionic and bosonic spins are contracted 
separately 



(El) 
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One might describe this vertex by the exchange of a real two-component boson field a = {(J^^a^) which couples to 
the composite operators ip''ip and ip'^ip with some Yukawa-type interactions. More explicit, we could use the following 
action: 

r.- I a.4,{Q)P„{Q)a^{-Q)- f Ki,{Q^,Q2)i^H-Qi)^{Q2)a^{Qi + Q2) 
JQ JQi,Q2 

haAQi^ Q2)^H-QiMQ2MQi + Q2). (E2) 



Qi 

Since u is a real boson, its propagator fulfills Pa-{Q) — P<j{—Q)- Together with Galilean invariance this implies that 
Pa does not have any frequency dependence Pcr{Q) = Pcr{ci) The exchange of a cr-boson corresponds to an 

instantaneous interaction. Note, however, that Galilean invariance is broken spontaneously by a condensate or a 
Fermi surface at nonzero density. In that case the cr-boson becomes dynamical and corresponds to a propagating 
phonon. An expectation value of a corresponds to a shift in the effective chemical potential 31]. 

In Eq. (|E2p the field a appears quadratic and we can eliminate it by solving its field equation. This results in: 

We observe that we get a tree- level contributions that has the spin-structure of the term ^ Asb. Assuming Yukawa 
couplings that are independent of the momenta and using the conventions of Sect. |V]the contribution to X^b reads 

—exchange 

^<7[Pl - P2) 



For an inverse propagator of the form P^iq) = rn^ + we find after the s-wave projection 

4piP2 \pI+pI-2piP2 



A3&,<T-oxchangc(Pl,P2;-E') = ^ m ( ^ -^^^ ] . (E4) 



The parameters and ha-^p , hcr^p can be chosen such that the form of Asf, is resembled closely. One can compare this 
to the tree-level contribution to Asf, by the exchange of a fcrmion ip. It is obtained from Eq. ([5]) by solving the field 
equation for and while using the fact that the propagator for -0 and the Yukawa coupling h are not renormalizcd: 



2h^ 

Pl^ + P2^ + (Pl + P2)^ - E- fl^' 



A3b,^-cxchangc(Pl, P2; -5') ~ „ 2 i „ 2 i i'^ , ^ n,2 'E' ~- (^^) 



After s-wave projection, this reads: 



\ ( p\ Jf_^ f pl+pl+PiP2-{t^^ + E)/2 \ 

A3&,V-cxchango(,Pl,P2;^j - -J^ m ^— g 7 ' ^^'^ 

2pip2 \pi +P2~ P1P2 - +E)/2J 

One can see, that the functional form of the two tree-level contributions after s-wave projection is quite similar. This 
is also the form of the momentum dependence found in the numerical solution of the flow equation for Aaf, without 
the cr-boson (Sect. I VIII and IVIII[) . We therefore expect that the description of Asf, as the exchange of a cr-boson with 
momentum-independent Yukawa couplings ha^p and h^^ leads to results that are comparable to the inclusion of the 
full (s-wave projected) momentum dependence. However, this description would be much more efficient with respect 
to the numerical effort. For the translation between the description used in the main part of this paper, where Aaf, is 
included as an own vertex and the description of Asb in terms of the exchange of a cr-boson, one might use the method 
of rebosonization [31 . 

The vertex Asa in Eq. ([5]) can also be described by the exchange of some particle, which corresponds in this case 
to a bound state of three atoms, the trimer or trion (2l| . Although the vertex Asf, and the momentum dependence 
of the Yukawa- like couplings were neglected in [2l|, the behavior found there was already qualitatively correct. Why 
this triatom approximation is not sufficient to describe the complete momentum dependence of A3 is discussed in 
Appendix [Cl 

We conclude that an effective description of the three-body physics with only a few couplings seems possible and 
would facilitate the study of systems at nonzero density and temperature. 
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